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ABSTRACT 

We present a new, BVIc photometric survey of the young open cluster IC 4665, which improves 
on previous studies of this young cluster by incorporating a rigorous standardization procedure, thus 
providing high-fidelity colors and magnitudes for cluster members. We use this new photometric 
dataset to reevaluate the properties (age and distance) of IC 4665. Namely, using a statistical approach 
incorporating r^ CMD modeling, we measure a pre-main-sequence isochrone age and distance of 36±9 
Myr and 360±12 pc, as well as a upper-main-sequence turn-off age and distance of 42±12 Myr and 
357±12 pc. These ages and distances are highly dependent on the isochrone model and color used for 
the fitting procedure, with a possible range of ~10-20 Myr in age and ~20 parsecs in distance. This 
spread in calculated ages and distances seen between colors and models is likely due to limitations in 
the individual membership catalogs and/or systematic differences in the predicted stellar parameters 
from the different sets of models. Interestingly, when we compare the isochrone ages for IC 4665 to 
the published lithium depletion boundary age, 28±5 Myr, we observe that this cluster does not appear 
to follow the trend of isochrone ages being 1.5 times smaller than lithium depletion boundary ages. 
In addition, comparing the overall magnetic activity (X-ray and Ha emission) in IC 4665 with other 
well studied open clusters, we find the observed activity distributions for this young cluster are best 
characterized by assuming an age of 30-40 Myr, thus in agreement with our pre-main-sequence and 
turn-off iaochione ages for IC 4665. Overall, although some age discrepancies do exist, particularly 
in the ages measured from pre-main-sequence isochrones, the range of possible IC 4665 ages derived 
from the various dating techniques employed here is relatively small compared to that found for other 
well studied open clusters. 

Subject headings: open clusters and associations: general — open clusters and associations: individual 
(IC 4665) — stars: evolution — stars: fundamental parameters 
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1. INTRODUCTION 

Open clusters have long provided necessary empirical 
constraints for our current models of stellar evolution 
due to their large collection of coeval, equidistant, and 
chemically homogeneous stars. Through observations of 
large numbers of open clusters with known ages, we are 
able to investigate the manner with which fundamental 
stellar properties evolve as stellar evolution progresses. 
This is particularly true for early stellar evolution (stel- 
lar ages < 100 Myr), when fundamental physical pa- 
rameters of stars change rapidly (e.g., radius, temper- 
ature, rotation rate, magnetic field strength). However, 
for these clusters to be used as valuable fiducial bench- 
marks of stellar evolution, we must be able to date them 
in a consistent manner with high-fidelity and precision. 
Currently, there are several dating methods utilized to 
determine open cluster ages, including fitting theoretical 
model isochrones to empirical H-R diagrams, light ele- 
ment abundance distributions, and exploiting the con- 
nection between magnetic activity and age. 

A heavily relied upon approach to determining 
absolute ages of star clusters is finding a best- 
fitting color-magnitude diagr am (CMD ) isochrone 
for th e nuclear turn-o ff (San daed 119581 : lMermilli"odl 
Il981al|bl : iMevnet et all ll993l ) and/or the pre-main- 
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sequence (PMS) tu r n-on (iD'Antona fc Mazzitellil 119971: 
iBaraffeet all 119981: iSiess et all I2000D to the cluster's 
main-sequence, typically performed using an ad-hoc 'by- 
eye' comparison of the model isochrone to the observa- 
tional photometric data. This technique, however eas- 
ily applied in theory, is in practice commonly ham- 
pered by uncertainties including field star contamina- 
tion, cluster distance and/or photometric error. Further- 
more, this technique suffers from the measured cluster 
ages being heavily dependent on the input physics used 
by the stellar evolution models. For instance, one can 
show for the same cluster that different nuclear turn- 
off and PMS mo dels can give stellar ages that differ by 
a factor of 1.5-2 (IChiosi et aLlHOOl IMevnet et al.l | 1993 : 



IMevnet fc Maederill997'. '200(i IBaraffe etal.lll998L 12002: 
Siess et al. 2000; Naylor 200§. 

There now exist several alternative dating techniques 
that employ other observable properties of stars to find 
absolute and relative ages for open clusters. Observa- 
tions have shown that the level of magnetic activity, in- 
dicated by the presence of X-ray and/or chromospheric 
activity, in solar-type stars clearly decreases wi th in- 
creasing age (Wilson 196 3; Skumanich 1972; Noves et al. 
1984; Sodcrblom et al. J1991I: iMamaiek fc Hillenbrand 
2008; .Cargile et al.. ,200i). This relationship is funda- 



mentally tied to the interplay between magnetic activ- 
ity and the rotationally-induced, magneto-hydrodynamic 
dynamo in solar-type stars. This causal relationship, 
the so called age-rotation-activity paradigm, generally 
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holds true; however, it can be affected by natural activ- 
ity cycles, similar to th ose seen in the Sun. For instance, 
iSoderblom et alj ()199lD show that uncertainties in stellar 
chromospheric ages can be as large as 50% for individ- 
ual stars. Of course, in open clusters this error can be 
statistically reduced by measuring many stars. 

IC 4665 is currently one of only a handful of well- 
studied PMS open clu sters (i.e., age of 5- 50 Myr), e.g., 
NGC 216 9 QO Mvr IJeffries et all [200l . NGC 2547 
(35 Myr Jeffries & Oliveiral 12005"), and IC 2391 (50 
Myr [Barrado y Navascucs et al. 2004). IC 4665 thus 
offers an interesting, additional opportunity to place 
empirical constraints on models of stellar evolution at 
young ages. Because of this, IC 4665 has been sur- 
veyed extensively to determine cluster membership, as 
well as observe physical properties of cluster mem- 
bers, including studies on the kinematics of the cluster 
(IProsseri 119931: IProsser fc Giampapa 1991 IManzi et al] 
120081: IJeffries et al.l I2009al). a ROSA T survev of stel- 
lar X-ray emission (jGiampapa et al.l I1998D . and stud- 
ies of the distribution of Li abunda.nces i n the cluster 
(IMartin fc Montesl[l997l: IManzi et alll2008l: [Jeffries et al.l 
l2009a| ). Most of these earlier investigations into the prop- 
erties of IC 4665 members relied on the o ptical photom- 
etry and membership catalogs reported in lProsseil ()1993l 
P93). Unfortunately, the cluster members detailed in 
this study rely upon uncertain photometric data, due to 
the use of an improper photometric standardization pro- 
cedure (see Section 12. 3|) , as well as having a high level of 
field star contamination as a result of their membership 
catalog being based upon unreliable proper motions (see 
Section IXTt. 

In light of such uncertainties, we report a new, stan- 
dardized BVIc photometric survey of the central 1 
sq. deg. of the cluster (see Section 12. 2|) . We use this 
new catalog to re-derive the age of IC 4665 using several 
different methods, including modeling the PMS (Section 
13. ip and upper-main-sequence (Section 13. 2|) photometric 
membership, and comparing the X-ray emission (Section 
13. 4|) and Ha emission fSection l3.5|) distributions to other 
well known open clusters. In addition, we compare our 
new age measurements in the context of the available 
ages for IC 4665 in the literature. 

2. THE OPEN CLUSTER IC 4665 

IC 4665 is a relatively young (< 100 Myr), nearby 
(^350 pc) open cluster in the constellation Ophi- 

uchus (R.A.(2000): 17''46™, Dec(2000):-h05°43'). The 
WEBDA Open Cluster Database0 reports for the clus- 
ter an age of 43 Myr, a distance of 365 pc, a high 
Galactic latitude {b = +17 deg), and has a reddening 
of E(B — V) — 0.174 (values originally published in 
iDias et aD 120021 ) . The youth and close proximity of this 
cluster has made it a popular target for studying stars 
during early times, while the majority of its low-mass 
members are still in the PMS phase. 

2.1. Extant Observations of IC ^665 

Previously, the most extensive op tical survey of 
IC 4665 was performed bv lProssen (|1993l ) , which included 

^ The WEBDA database, developed by J.-C. Mermilliod, can bo 
found at http://www.umvie.ac.at/webda/ 



astrometric, photometric, and spectroscopic programs to 
identify members of the cluster. P93 acquired BVRIfc 
CCD photometry for 4100 stars down to a limiting mag- 
nitude of T^ '^ 18.5. Proper motions were measured from 
archival photographic plates, and membership probabil- 
ities were derived for stars brighter than V~14.5. Their 
analysis found 170 stars with non-zero proper motion 
membership probability, of which 22 were identified as 
previously known cluster members, 57 were classified as 
photometric candidates {i.e., stars lying near a 36 Myr 
isochrone), 72 photometric nonmembers, and 22 stars left 
with undecided membership. They also obtained low- 
dispersion spectra for 41 faint candidate stars, of which 
25 showed definite Ha emission while 5 showed possible 
weak Ha emission. P93 include these 30 stars as proba- 
ble members of IC 4665. 

Multiple investigations into chromospheric and coro- 
nal activity have been conducted in IC 4665. Ha emis- 
sion, as an indicator of chromospheric activity, has also 
been observed in 14 F-^ G-, and K-spectral type stars 
(jMartin fc MontesI 119971) . as well as for 3 9 lower-mass, 
M-spectral type stars (P93). More recentlv. IJeffries et al.l 
([^009a) measured Ha equivalent widths for 56 IC 4665 
stars ranging from early-F to mid-M spectral type. In 
addition to chromospheric studies, IC 4665 also has been 
observed with the ROSAT X-r ay telescope, in o r der to 
detect coronal X-ray emission. ICiampapa et al.l (|1998[ ) 
reported the results of a 76.9 ks exposure of the cluster 
with the HRI instrument on ROSAT, from which they 
found 28 X-ray detections matching optical counterparts 
in the P93 catalog. Six (6) of these 28 X -ray sources were 
shown to be cluster nonmembers by iGiampapa et all 
()1998() based on spectral-type classifications. 

IC 4665 is currently one of just five open clusters 
with a detected lithium d epletion boundary (hereafter 
LDB). IManzi et all ()2008f ) measured an age of 28 ± 5 
Myr for the cluster based on the absolute magnitude of 
stars found at the LDB. Moreover, the global lithium 
abundance distributio n in I C 4665 stars was first studied 
by Martin & Monteg (jl997t ). where they measure the Li 
abundance of 14 early-G through early-M spectral type 
stars. Interestingly, they found that the Li distribution 
of IC 4665 is similar to the older Pleiades open cluster 
(~125 Myr), although they had a limited sample size of 
measured Li abundances. More recently, IJeffries et al.l 
(|2009a|) determined the Li abundance of 20 late-F to 
early-M spectral type stars, and confirmed that the Li 
distribution of IC 4665 appears more like the Pleiades 
than other younger open clusters. 

The stellar rotation distribution of IC 4665 is com- 
monly used to empirically constr ain theories o f early 
angular momentum evolution (e.g. iBarnesI [20031 [2007). 
Using rotation period measure ments for nin e solar-type 
stars (F9 to KO spectral type), lAUain et all fl996 ) find 
that there is a dearth of slow rotators (periods > 4 days) 
in IC 4665, indicative of a young cluster w here the stars 
have no t had time to spin down. However. IJeffries et al.l 
((2009a) measure projected rotation velocity (Vsinz) for 
IC 4665 cluster members, and find that there is a lack of 
fast rotators (Vsini > 20 km s^^) with 1 < (F - /c) < 2 
(~ K-type stars). Comparing this Vsini distribution to 
the older IC 2391, IC 2602, and Pleiades, they find that 
the K-type stars in IC 4665 appear to rotate more slowly 
than one would expect at ~30 Myr. Finally, a recent 
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survey reported rotation per iods for 20 IC 466 5 mem- 
bers with masses < 0.5Mq (jScholz et al.ll2009D . They 
show that all of these low-mass stars are rotating rapidly 
(period < 1 day) , a finding similar to studies of other 
PMS clusters fe.g. llrwin etalll2008n . 

2.2. New SMARTS BVIc Observations 

We observed IC 4665 using the 1.0m SMARTS tele- 
scope (the old YALO telescope) situated at the Cerro 
Tololo Interamerican Observatory (CTIO), Chile, during 
the nights of 15, 16, 17 September, 2005. The data were 
acquired with the quad-amplifier Y4KCam CCD camera, 
equipped with Johnson-Cousins BVIc filters, with its 15 
/im pixels and 0.289 arcsec/pixel plate scale results in 
an on-sky field of view per CCD field of 19f3 x 19.'3. 
Nine fields were observed covering a uniform grid of a 
3x3 mosaic of CCD fields covering approximately 1 
square degree, centered on R.A.(2000):17M6™13.5' and 
Dec(2000):-f05°41'52". We employed exposure times of 
100, 50, and 25 seconds for the B-, V-, and Ic-band ob- 
servations, respectively. The central coordinates of each 
IC 4665 field for which BVIc photometry was obtained 
are shown in Fig. [Hand summarized in Tabled] 



TABLE 1 
Summary of SMARTS 1.0m Observations of IC 4665 
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Fig. 1. — Digital Sky Survey image of the areal region covered by 
our SMARTS IC 4665 survey. Each square represents an observed 
field (coordinates given in Table [TJ. Also plotted is the overlapping 
control field (CF) identified by the red dashed bo x, and the cluster's 
core radius (black dashcd-circle) as derived by IKharchenko et aP 
||200S '|. 



All science and standard star images were processed 
for overscan region subtraction, master-bias subtraction 
and fiat fielding, using twilight sky images, using stan- 
dard procedures in the IRAlfl suite of data reduction 
algorithms. Substantial use was also made of two batch 
processing IRAF scripts, written by Phil Massey of Low- 
ell Observatory, whose protocol specifically handles the 
FITS headers and quad-amplifier readout of the CCD 
camera. The relatively low space density of the cluster, 
and its lack of substantial field crowding, allows us to 
employ aperture photometry for both IC 4665 science 



IC 4665 


R.A.'' 


Dec^ 


Date of 


Field^ 


[HH:MM:SS] 


[DD:MM:SS] 


Observation 


Fl 


17:47:29.8 


-(-06:01:00.0 


16/Sept/2005 


F2 


17:46:13.5 


-1-06:01:00.0 


16/Sept/2005 


F3 


17:44:56.3 


-1-06:01:00.0 


17/Sept/2005 


F4 


17:47:29.8 


-(-05:41:52.0 


15/Sept/2005 


F5 


17:46:13.5 


-(-05:41:52.0 


15/Sept/2005 


F6 


17:44:56.3 


-(-05:41:52.0 


17/Sept/2005 


F7 


17:47:29.8 


-(-05:22:56.0 


15/Sept/2005 


F8 


17:46:13.5 


-(-05:22:56.0 


16/Sept/2005 


F9 


17:44:56.3 


-(-05:22:56.0 


16/Sept/2005 


CF 


17:46:52.6 


-(-05:32:00.0 


15/Sept/2005 


CF 


17:46:52.6 


-(-05:32:00.0 


17/Sept/2005 



^ Coordinates are Y4KCani field centers. All coordinates 
are J2000.0 epoch. 

^ F'n'-fields represent the main survey area of the cluster 
in a 3 X 3 mosaic covering about one square degree centered 
close to the cluster's WEBDA coordinates. Photometric 
data for a control field (CF) were also procured to allow 
us to assess the photometric integrity of the dataset. 

and standard star images. Source searching and aper- 
ture photometry, with an aperture radius of 13-pixels, 
were achieved using the DAOPHOT package in IRAF 
IStctson 1987, 1992). 

Equatorial BVIc standard stars cataloged in iLandoltl 
(fT992.. .2009. ) , with a broad dynamic range of photomet- 
ric magnitude and color, and located at wide distribution 
of airmass, were observed nightly in order to correct for 
atmospheric extinction and to transform extinction cor- 
rected, instrumental magnitudes onto the standard sys- 
tem. Modified forms of Bouguer's law, see Eqs l2.2l fa-d). 
were used to define the relationship between standard 
and instrumental magnitudes for the standard stars. 



V = v + e{B~V)+Cv-ii^ X (1) 

V = v + e{V-Ic)+ Cu -nv X (2) 
{B-V)^ ixib -v)+^bv-fi Kbv X (3) 

{V - Q = ^{v - Zc) +U-^ i^vi X (4) 

where 

• V, {B — V) and {V — Ic) are magnitudes/indices on 
the BVIc standard system 

• v, (b — v), and {v — ic) are measured, instrumental 
magnitudes/indices 

• Ky,Kiiy and Kyi are filter dependent extinction co- 
efficients 

• e, fj, and V' a-re color transformation coefficients 

• Cf I Cfcu and £,m are zero-point coefficients. 



• X = airmass of target 
zenith distance 



sec z - where z is the 



** IRAF, in our case through http://iraf.net is distributed by the 
National Optical Astronomy Observatories, which are operated by 
the Association of Universities for Research in Astronomy, Inc., un- 
der cooperative agreement with the National Science Foundation. 



Unknown extinction coefficients, color transformation 
coefficients and zero points were determined by solv- 
ing the self-similar series of linear simultaneous alge- 
braic equations detailed in Eas l2.2l (a-d). Solutions were 
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TABLE 2 

Nightly Extinction, Transformation and Zero-point Coefficients 



Date 



Equation Extinction'' 
Coefficient 



CTC*^ 



Zero 
Point 



Mean'' 
A mag 



RMS'* 
A mag 



No. of= 
Standards 



15/Sept/2005 


V (B-V) 


0.1478 


0.0140 


-1.8654 


6.3 X 10-05 


0.0148 


90, 


72 




V (V-Ic) 


0.1560 


0.0128 


-1.8541 


1.0 X 10"°^ 


0.0176 


87, 


75 




B-V 


0.1116 


0.8940 


-0.3011 


4.5 X 10^05 


0.0166 


90, 


71 




V-Ic 


0.1235 


1.0009 


1.0472 


1.0 X 10^05 


0.0172 


89, 


58 


16/Sept/2005 


V (B-V) 


0.1476 


0.0166 


-1.8518 


2.4 X 10-05 


0.0157 


80, 


62 




V (V-Ic) 


0.1492 


0.0171 


-1.8524 


1.0 X 10-05 


0.0170 


74, 


62 




B-V 


0.1497 


0.8731 


-0.2334 


7.5 X 10-05 


0.0177 


80, 


60 




V-Ic 


0.0832 


1.0001 


0.9935 


6.5 X 10-05 


0.0172 


80, 


60 


17/Sept/2005 


V (B-V) 


0.1556 


0.0272 


-1.8307 


8.1 X 10-05 


0.0180 


76, 


57 




V (V-Ic) 


0.1574 


0.0239 


-1.8290 


-5.0 X 10-05 


0.0183 


76, 


57 




B-V 


0.0460 


0.8808 


-0.3565 


-7.7 X 10-05 


0.0220 


64, 


52 




V-I, 


0.0163 


1.0040 


0.9107 


3.3 X 10-05 


0.0190 


74, 


54 



• ■ '■■■ ';S^:is--^'>w:^;,^;--:s.. 



^ In units of magnitude/airniass. 

^ Color Transformation Coefficient 

'^ Observed mean difference between calculated and published magnitudes or colors for ILandolj 

{1993, 120091 ') standard stars. 

* * O bserved RMS difference between calculated and published magnitudes or colors for ILandolj 

||M)_2, 2009) standard stars. 

^ m, n - where m is the number of Landolt standard stars observed, and n represents the number 

of standard stars used in the final photometric solution. 

derived using a least-squares fit algorithrr0 which pro- 
cesses the numeric parameters via an iterative refinement 
method. Comparing calculated magnitudes and colors 
of Landolt standard stars to their published values al- 
lows an iterative rejection process to take place, thereby 
eliminating seriously discrepant standard star measure- 
ments - of which there were few. Solutions to the color 
equations for each observing night are listed in Table [21 
whereby these solutions resulted in differences between 
calculated and published magnitudes of Landolt stan- 
dard stars that were never greater than 2.2% {i.e., 0.022 
magnitudes). Target star photometry was readily deter- 
mined by substituting the derived extinction, transfor- 
mation and zero-point coefficients along with measured 
instrumental magnitudes into Eas l2.2l fa-d). 

The full BVIc catalog for all IC 4665 fields is pro- 
vided in Table [S] including their associated 2MASS J, 
H, and Kj photometric and astrometric data, corre- 
lated with point sourc es in the online 2MASS database 
(jSkrutskie et al.ll2006l ). Of the 7359 stars in our mas- 
ter BVIc catalog, we find 452 stellar sources with- 
out matches within 3.0 arcseconds of 2MASS sources. 
For these stars, we provide J2000 astrometry, precise 
to ±1 arcsecond, using a 6th-order polynomial fit to 
X,Y CCD coordinates and accurate and precise coor- 
dina tes for reference sta rs in the SuperCOSMOS cata- 
log (|Hamblv et al.|[200lD . A visual presentation of the 
BVIc photometric catalog is shown in Fig. [21 in which 
we p lot V versus B — V and V — Ic CMDs in concert 
with D'Antona fc Mazzitelhl (|1997f ) fiducial isochrones. 

In order to investigate the stability of our photometric 
system, we also observed a control field (CF) in IC 4665 
during the nights of the 15*^ & 17*^ of September 2005. 
Fig.[3]shows our V-, B-, and Ic-band magnitude compar- 
isons for stars in common to the control fields observed 
on the two different nights. A simple statistical analy- 
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Fig. 
metric 



5 The algorithm, F04AMF, 
gorithms Group (NAG) 



is distributed by the Numerical Al- 



2. — Color-magnitude diagrams for our full photo- 

catalog ar c plotted (black dots). Also provided are 

ID'Antona &: MazzitcUi (1997) isochrones (blue dotted line) for the 
cluster's lithium depletion boundary age (27. 8 Myr; [M anzi ct aT] 
l200ai') and HIPPARCOS distance (355 parsecs; Ivan Leeuweni2009i r 

sis of these data is presented in the Table [H From these 
plots and table, we can infer that the V- and Ic-band zero 
points are at least consistent at the level of ^ 3 — 4% for 
objects with V and Ic <16 magnitudes. For these same 
objects, the B-band observations do however show an 
offset at the ^ 5 — 6% level. Naturally, for the dimmer 
objects (V, B, Ic >16 magnitudes), with decreasing S/N 
levels, diminished precisions of order '--^ 9 — 10% are seen. 

2.3. Comparison to Prosser (1993) Photometry 

Existing studies reporting on the physical character- 
istics of IC 4665 solar-type stars are under-pinned on 
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SMARTS BVIc AND JHKs 


Photometry 


FOR THE Region Around IC 4665 






Mosaic Object 
Field Number 


R.A. Dec. 

[HH:MM:SS] [DD:MM:SS] 


V 


B-V 


V-lc 






F6 
F6 
F6 


682 
259 
723 


17:44:17.47 
17:44:17.59 
17:44:17.62 


+05:38:34.2 
+05:40:44.7 
+05:36:05.0 


18.299+0.193 
15.961+0.022 
17.051+0.061 


99.999+9.999 
1.083+0.063 
99.999+9.999 


1.566+0.243 
1.232+0.033 
1.006+0.101 












TABLE 3 

CONT. 








Mosaic 
Field 


Object 
Number 


R.A.(2MASS) 
[HH:MM:SS] 


Dec.(2MASS) 
[DD:MM:SS] 


Pes. Off. '^ 


J 


H 


K, 


2MASS 
Qual Flag 


F6 
F6 
F6 


682 
259 
723 


17:44:17.41 
17:44:17.55 
17:44:17.59 


+05:38:34.60 
+05:40:44.52 
+05:36:04.49 


0.96 
0.60 
0.62 


16.149+0.108 
13.965+0.030 
15.437+0.071 


15.819+0.167 
13.457+0.036 
15.027+0.088 


15.470+0.209 
13.325+0.040 
15.181+0.161 


ACC 
AAA 
AAC 



Note. — Only a few rows are displayed here to show the content of the table, the full table is available in electronic form. All 
positions are given in J2000 coordinates. Magnitude, colors, positions, and uncertainties given as 99.999 or 9.999 indicate that 
the value is undetermined. 

'^ The positional offset between the R.A. and Dec. given in our photometric catalog and the corresponding 2MASS source 
position. 



TABLE 4 
Internal Precision of SMARTS BVIc Photometry 



Filter 


15th - 17th 
Mean Num 


15th - 16th 
Mean Num 


16th - 17th 
Mean Num 


Full Catalog 


V 
B 

Ic 


0.024+0.056 
0.079+0.061 
0.029+0.071 


1106 
669 
1058 


0.019+0.062 
0.058+0.070 
0.064+0.061 


234 
107 
219 


0.030+0.066 
0.032+0.075 
0.004+0.076 


206 
110 
190 


V < 16 


V 
B 

Ic 


0.016+0.014 
0.066+0.038 
0.020+0.014 


365 
166 

794 


0.030+0.014 
0.054+0.040 
0.044+0.013 


72 
31 
155 


0.001+0.015 
0.020+0.045 
0.011+0.015 


65 
30 
138 


V > 16 


V 
B 

Ic 


0.028+0.077 
0.092+0.124 
0.036+0.091 


741 
503 
264 
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0.010+0.098 


141 
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the P93 photometry study, and it is therefore imperative 
that we determine the rehabihty of their photometric sys- 
tem as compared to ours. A detailed examination of the 
standardization procedure of P93 shows that, although 
nightly standard stars were used to derive extinction co- 
efficients and color-transformation equations, the zero- 
points of their BVIfc photometric system are tied to a 
solitary blue, high-mass B9e star. This star, HD 161261, 
has subsequently been fou nd to have V-band variabi l- 
ity on the order of 0.1 mag (jKazarovets fc Samuslll997D. 
which maybe due to a binary component (|Pourbaix et all 
l2004f l. as well as an infra.red e xcess typical of rapidly 
rotating Be stars (jYudinI 1200 It ) . Not surprisingly, we 
observe systematic offsets between the P93 photometry 
and our own that are likely a result of this extrapolation 
of standardized photometry derived from such a single, 
blue, possibly variable and binary standard star. 

For stars in common to both catalogs, we plot in Fig.|4] 
comparisons between the P93 study and our own V, 
B — V, and V — Ic photometric datasets. We trans- 
lated the P93 Kron V — I colors to the Cousins system 



using the formalism given in iBessell fc WeisI (119871 ) . We 
see a significant amount of scatter in the P93 photome- 
try (^0.1 mag in V; ^0.2 in i? — V/V — Ic) compared to 
our optical catalog. In addition, there appears to be a 
systematic offset of '--^0.05 mag in V, B — V, and V — Ic 
in the faintest /reddest P93 objects. Such a magnitude of 
photometric scatter and systematic offset may cause sig- 
nificant problems when trying to calibrate both color and 
magnitude relationships for stars in IC 4665; therefore, 
for the purpose of this manuscript, we will preferentially 
employ our, new BVIc photometry in all subsequent anal- 
ysis. 

3. THE AGE OF IC 4665 

Determining a precise and accurate age for IC 4665 is 
fundamental to its use as a constraint, or boundary con- 
dition, for theoretical models describing early stellar evo- 
lution. Its age however, continues to be a matter of un- 
certainty. Ages for IC 4665 found in the literature range 
from greater than 25 to less than 100 Myr, although one 
of the age determinations is derived using the lithium 
depletion boundary method (|Manzi et al..,2008.) . which 
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Fig. 3. — Comparison of CTIO BVIc photometry for control 
fields in IC 4665, based on observing nights 15, 16, and 17 Sept. 
2005. Upper panels show direct comparison, while their residuals 
are plotted in each lower one. Data arc plotted with their associ- 
ated photometric error. Red dotted lines are loci of equality and 
not fits to the data. 

SO far has only been applied to a small handful of open 
clusters. With such a large range of possible ages cited 
for IC 4665, its applicability for constraining theoreti- 
cal stellar evolution models remains suspect. Here, we 
use our new photometric data for cluster members of IC 
4665 to take a new look at its discrepant age determina- 
tions to verify whether there exists a true, observational- 
dependent age spread in IC 4665, or whether this spread 
in ages is due to uncertainty in the extant studies of the 
clusters. 

3.1. Pre-Main-Sequence Age 

In the P93 study, the authors reported that they were 
unable to determine the age of IC 4665 from the location 
of the low-mass pre-main-sequence, due to the lack of a 
defined low-mass sequence in their CMD. This may have 
been the result of systematic scatter in their photome- 
try and/or their cluster membership determination being 
based on unreliable proper motion data. More r e cently , 
the infrared study of IC 4665 by Ide Wit et all (pOOQ ) 
suggests that the cluster has an age of 50 to 100 Myr 
based on the location of its low-mass cluster sequence. 
However, Ide Wit et al.l note that an accurate age for the 
cluster is difficult to measure due to the large amount of 
field star contamination present in their survey. 

A quick examination of the full IC 4665 CMD reveals 
that we too confirm that an accurate fit of the cluster 
sequence is particularly difficult to achieve due to field 
star contamination in our color-magnitude diagram (see 
Fig. [2]) . While IC 4665 lies relatively far from the Galac- 
tic plane {b = -1-17°) the level of contamination along 
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Fig. 4. — V, B — V, and V — Ic residual comparison plots for 
the P93 catalog and our optical catalog (CJ). Our translation from 
P93 V - Ik to V - Ic IS described in Section [23] The blue dotted 
lines indicate equality between the systems, whereas the red dashed 
lines are calculated means for V < 14.5 & > 14.5; B — V < 1 &c > 
1; V — Ic < 1.5 & > 1.5, respectively. 

the CMD suggests that there is a high stellar density 
in the cluster's line of sight through the Galaxy. In or- 
der to account for the field star contamination in the 
IC 4665 CMDs, we further constrain our analysis to those 
stars having ancillary membership indicators from a dis- 
parate range of extant photometric and spectroscopic ob- 
serving campaigns. In the first instance, we position- 
ally matched our photometric catalog with data from 
the following sources: three catalogs of stars with com- 
mon kinemati c properties (proper rnotions -- P93: radial 
velocities — iProsser fc Giampapal 119941: iJeffries et alj 
l2009al) . stars with X-rav emission iGiampapa et al.lll998l 
two catalog s of stars with significant lithium absorption 
(jMartin fc M ontes 1997; Jcffi-ies ct al. 2009j.), three cat- 
alogs of stars exhibiting strong Ha emission (EW Ha > 
Ol — P93: lMartin fc MMrtesl[l997l: IJeffries et al.l[2559al) . 
and finally, two catalogs of stars with short rotation 
periods in the re gion of IC 4665 (jAllain et al.l 119961 : 
iScholz et al.ll2009l ). X-ray and Ha emission (coronal and 
chromospheric) , the presence of lithium in stellar pho- 
tospheres, and rapid rotation are all indicators of youth 
for low-mass stars, and therefore we can exploit such sur- 
veys in order to segregate young cluster members from 
the typically much older field star populations. In Fig.jS] 
we plot the CMDs with stars identified from the vari- 
ous membership catalogs. We will subsequently exclude 
proper-motion-selected members from our cluster cata- 
log due to the lack of a clear (pre-) main-sequence in its 
CMD (see Fig. [5] top panels), viz, proper motion must be 
an extremely poor membership criterion for this cluster. 
This can be explained by the similarity of the cluster's 
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systemic proper motion to that of the general field pop- 
ulation, justifying the exclusion of the cluster's proper 
mot ion catalog from all subsequent analysis of IC 4665 
(see lJeffrieset al. 2009a). 

We recognize that the membership catalogs cited above 
contain varying levels of field contamination, as well as 
biases based on their color and/or magnitude ranges. For 
instance, although our catalog of radial velocity mem- 
bers contains the largest number of stars spread over 
the largest range of colors, it also has a high probabil- 
ity of field star contamination due to the similarity of 
the kinematic distribu tion of the clus ter compared to 
that of the field (.Jeffries et al.ll2009aD . On the other 
hand, the catalog of stars showing X-ray emission has 
a much lower probability of contamination due to the 
ubiquitous nature of high magnetic activity in young 
stars compared to old field stars. However, the rela- 
tively low number of stars identified can be attributed 
to the fairly narrow field-of-view of the ROSAT X-ray 
satellite compared to the on-sky distribution of appar- 
ent cluster members. In an attempt to account for the 
relative po wer of using multiple membership criteria, 
iJeffries et al. (2009a) produced a catalog of highly proba- 
ble members of IC 4665 using a combination of kinematic 
motions, lithium content, optical and 2MASS photom- 
etry, and spectral-index indicators as filters to remove 
foreground and background field star interlopers in the 
cluster catalog. We include this list of high-fidelity clus- 
ter members in our analysis, and note that it is most cur- 
rent, up-to-date compilation of bona fide IC 4665 mem- 
bers. 

In an effort to find an initial estimate for the age of 
IC 4665 using its low-mass population, we plot in Fig. |6] 
the V versus V — Ic color-magnitude diagram for con- 
firmed members of three open clusters; for IC 4 665, the 
membership catalog from IJeffries et all (|2009al ) is used, 
whereas the 35 Myr old NGC 2547 cluste r, we employ 
the ph otometric membership catalog from iNavlor et aLl 
(|2002D . Finally for the -100 Myr old Blanco 1 cluster, 
James et al. (m prep) provide us with its photometric 
catalog of kinematic members. While some photomet- 
ric scatter is apparent, probably due to binarity and 
a small number of cluster non-members, these CMDs 
evidently show that the lower mass stars in IC 4665, 
{V — Ic)o > 1.00, lie systematically above the Blanco 1 
main sequence. In fact, in this color range IC 4665 
stars have colors and magnitudes closely matching the 
corresponding stars in NGC 2547, suggesting that at 
least photometrically, IC 4665 is of a comparable age to 
NGC 2547, i.e., ~ 35 Myr. Although this age estimate 
is insightful, we now establish a more robust determina- 
tion of the cluste r's age using the so-call ed r^ isochrone 
fitting technique (jNavlor fc Jeffriesll2006( ). 



3.1.1. 



Isochrone Modeling 



The classic approach in determining fundamental prop- 
erties of a young open cluster, such as age and distance, 
is finding a best-fitting CMD isochrone for the apparent 
cluster sequence. Typically, this is performed using a 
"by-eye" comparison of the model isochrone to the ob- 
servational data, which can be burdened by uncertainties 
due to field star contamination and/or photometric er- 
ror, as well as an apparent spread in the cluster sequence 
due to binarity. In this manuscript, we describe a new. 



statistical approach in modeling the lower-main-sequence 
of IC 4665. 

The Tau-squared (r^) maximum likelihood statis- 
tic has recently been developed for modeling 
CMDs with two - dimen s ional t h eoret ical isochrones 
(jNavlor fc JeffriesI [20061 : iNavloil [20091 ). Examples 
of its application can b e found in iJefjF ries et ^ 
(|2007[): iJoshi et al.l (|2008[) : iMavne fc Navlori (,200^ 



IJeffries et al.l ()2009b[ ). In short, t^ is a generalized x^ 
statistic that takes into account correlated uncertainties 
in two dimensions. Thus, we can use this statistical 
tool to effectively model a two-dimensional distribution, 
in our case an age- and distance-dependent CMD 
isochrone, to the V/B ~ V and V/V - Ic CMDs by 
minimizing r^ to find a best-fitting isochrone. 

Following clo s ely th e technique described in 
INavlor fc JeffiTe^ (|2006l ): iNavlori (pOOl ). we produce a 
grid of stellar model isochrones spanning a range of ages 
and distances. We use the low -mass, solar-metallicity 
stellar evolution models fro m iD'A ntona fc Mazzitellil 
(fT99l . iBaraffe et all (fl998l ). and~[Sicss et al. (2000), 
and employ a Pleiades tuning method to translate 
model luminosities and effective temperatures (T^ff) 
into magnitudes and colors (for a full descript i on of 
the Pleiades tuning method, see IJeffries et al.l 120011 : 
INavlor et aLll2002D . The r^ softwareO uses Monte Carlo 
methods to produce a probability, based on r^, of 
finding a star drawn from the full photometric catalog 
(including photometric binaries) represented by each 
age-distance pixel on the grid. These probabilities are 
calculated for all stars in the CMD, and summed to 
derive the overall r^ for the given age and distance 
isochrone. By repeating this process for a range of ages 
and distances, a r^ grid is produced, and the minimum 
value is the best-fit isochrone for the modeled CMD. 

In order to have our best-fit isochrones not be influ- 
enced by obvious outliers, we added an additional soft- 
clipping rejection step (Na vlor fc Jef fries 2006), where 
by points having colors/magnitudes lying several a away 
from the range of model isochrones are ab initio removed 
from the modeling process. In our modeling of IC 4665, 
we assume an extinction relation of A„ = 3.12xE{B—V), 
a reddening of E{B -V) = 0.174 (WEBDA Open Clus- 
ter Database), an E{V - Ic) = 1-25 x E(B - V) rel ation 
and a binary fraction of 0.5 (jCrampton et al.lll976( ). 

We have assumed in our analysis that differential red- 
dening is neg ligible for IC 4665 . Th is is based on 
the results of [Crawford fc BarnesI (|1972D study finding 
the reddening in B — V varied by only —0.02 magni- 
tudes for 17 high-mass (B spectral- type) stars across the 
cluster. This amount of variation in B — V color ex- 
cess is at the level of our photometric internal precision 
(A(B - y) - 0.017 - 0.022 mag) as determined by our 
nightly comparison analysis (see Section 12.21 and Table 

Once a best-fit isochrone for each of our IC 4665 CMDs 
is found, we derive the probability of the isochrone model 
being a good fit based on their minimum r^ values. We 
include an additional magnitude- independent systematic 
uncertainty to the colors and magnitudes in our CMDs 
to increase their probabilities to an acceptable value 

^ Th e T^ software is publicly distribut ed by T . Naylor at 
|http://www.astro.ex.ac.uk/pcoplc/tirrm/tau-squared[ 
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Fig. 5. — Intrinsic V — Ic and B — V color-magnitude diagrams for IC 4665 are presented (object density given by gray contours) with 
cluster cand idate members identified (crosses) using a secondary membership criterion, labeled in the lower-left corner of each panel (see 
Section [XT] for source references). We use a cluster reddening of E(B — V) = 0.174 (WEBDA), an E(y — I) = 1.25 X E{B — V), and an 
Ai, = 3.12 X E{B — V) relation, in order to derive intrinsic colors and magnitudes. Colors in each panel correspond to the associated data 
points in Figure |9] 
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Fig. 6. — Absolute V versus intrinsic V — Ic magnitudes are 
plotted for the open clusters IC 4665, NGC 2547 (35 Myr), and 
Blanco 1 (100 Myr). Intrinsic colors and magnitudes were cal- 
culated using cluster's E{B - V) data (WEBDA), an A„ = 
3.12 X E{B - V) and an E{V - h) = 1.25 X E{B - V) relation. 
The apparent-to-absolute magnitude conversion was achieved using 
the ap propriate HI PPA RCOS distances reported in van Lecu wer] 
1120091) . See Section lS.ll for data sources. 

(i.e.,> 30%). This is analogous to the traditional method 
of increasing error bars to arrive at a x^ ~ 1, and ensures 
that the uncertainties on our ages and distances are ap- 
propriate relative to the scatter in the IC 4665 CMDs. 
Typically, we found that an additional ^ 0.01 — 0.03 
increase is required for our V, B, and Ic magnitudes 
to achieve sufficiently high best-fit probabilities. This 
magnitude of systematic uncertainty can be easily ex- 
plained by sys tematics in the standardization procedure 
(jNavlor et al.1 12002') or natural variability in such young, 
solar- type stars {e.g., spotted stars in IC 4665 have re- 
ported optical va riability amplitudes of ~ 0.03 — 0.1 mag, 
lAllain et al.][l996l'l . 

As an example of the basic procedure behind r^ 
isochrone modeling, we plot in Fig. [7] the r ^ grid derived 
from V — Ic CMD isochrone modeling of the [Jeffries et aLl 
(l2009afl membership ca talog using the PMS models from 
iD'Antona fc Mazzitelli The displayed color-scale indi- 
cates the total T^ values for the indicated distance mod- 
ulus (hereafter DM) and age grid points. The minimum 
value in this grid defines the best-fit age and DM for this 
membership catalog. This grid of r^ values can also be 
thought of as a grid of probabilities (see above); there- 
fore, r^ values can be found for different confidence levels, 
e.g., displayed in Fig. [7] is the 68% confidence contour. 
The ranges of ages and distances falling within this con- 
tour defines the uncertainty in our best-fit measurements. 
We present in Fig. \8\ the best-fit isochrone from this t^ 
grid analysis. 

3.1.2. Pre- Main- Sequence Isochrone Results 

In Fig. [9] we show the best-fit ages and distances, de- 
rived from our Tau-squared fitting procedure, for each of 
our seven IC 4665 membership catalogs (5 individ ual cat- 
alogs described previously, the catalog from Jeffri es et alJ 
(|2C)09ai) ■ and a combined catalog that includes any ob- 
ject with some type of membership criterion), and three 
stellar evolution models representing two different colors 
{B — V and V — Ic). We note that we were unable to 
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Fig. 8. — C olor-magnitude diagram for the objects in 
the IJeffr ies et at] membership catalog. Plotted is the best-fit 
ID'Anton a & Mazzitelli isochrone from our r^ analysis, with the 
color-scale indicating the simulated relative number density of stars 
along the cluster sequence. The red circles indicate objects that 
were rejected by our "soft-clipping" scheme, and the blue circles 
are the data used in the isochrone modeling. 
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Fig. 9. — Best-fit ages and intrinsic distance moduli resulting from our r^ isochrone modeling of the low-mass members of IC 4665. 
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by the legend in the center of the plot (see Section |3. II for source references). Along the left of each plot, the final number of data points 
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on the left. For comparison purposes, we plot a black dashed-li ne box in each box, wh ich represents the approximate range of ages and 
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find a best-fit age and DM solution for every membership 
catalog using each model and/or color, with the rotation 
period and Hq catalogs being particularly problematic. 
The failure of these model convergences is primarily due 
to the specifics of each membership catalog {e.g., large 
amounts of field-star contamination, not having a wide 
color-magnitude range in data points, etc.). 

In the following analysis we choose to use the member- 
ship catalog from Jeffries et al. (2009a) as the exemplar 
best-fit age and DM for the different models and colors, 
as it has the lowest probability of contamination by field 
stars (see Section [37T]) . Nevertheless, modeling multiple 
membership catalogs allows us to identify the possible 
range of ages and DM one should expect when using spe- 
cific membership criteria to identify cluster members. In 
Tableini we list the median of the ages and distances, de- 
rived for all membership criterion catalogs, together with 
the observed range of values from our best-fit isochrone 
for each color and stellar model. We also include the 



best-fit solution for the I Jeffries et al.l membership list. 

For the most part, the best-fit ages and DM found for 
all three PMS models are typically within their individ- 
ual errors of each other, suggesting that statistically they 
are all three are in agreement of IC 4665 's age and dis- 
tance. However, we do find that some points lie beyond 
1 a from the other derived ages and DM (e.g., the best- 
fit para meters found for the V — Ic ^'i catalog using the 
iD'Antona &: Mazzitellil models), and are most likely the 
result of a poor fit due to limited or discrepant data. 
This is evidence for the strength of modeling multiple 
membership catalogs, that we are able to clearly adjudge 
whether an individual catalog's measured age and DM is 
discrepant due to a poor isochrone fit or otherwise. 

Here, we would like to highlight some results of our 
PMS isochrone modeling of IC 4665. First, if we focus 
solely on the B — V isochrones, we find good agreement 
in the measured DM for all three mode ls, w hereas only 
the ages for the ID'Antona fc Mazzitellil and iSiess et aD 
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TABLE 5 
IsocHRONE Ages and Distances for IC 4665 



Parameter'' 



Median'' 
B-V V - Ic 



Jeffries et al. 

B^v — r^7. 



DAMdr" 



Age [Myr] 
D.M. [mag] 
(Dist. [pc]) 



48.0±11.0 
7.72±0.08 
(350±13) 



38.5±16.5 
7.76±0.19 
(370;* 



^"--ID 



10 

1 

+0.03 
0.05 

(358l|) 



39.0 

7.77 



+3 

-1 
+0.04 
-0.07 

+ 7 



68. 7f^ 

,+0.05 
-0.01 

(350t«) 



32 

7.79 
(361+{i) 



-5 

7 oo+O.ll 
'•»»_0.01 
,+ 19) 



BCAH98'^ 



Age [Myr] 
D.M. [mag] 
(Dist. [pc]) 



65.2±14.1 

7.72±0.05 

(350±8) 



31.5±10.5 

7.94± 0.10 

(387t^?) 



7.72"* 



TTD- 



39. ( 



(377: 



27.5Zr 
7 06+003 

(39ll?i) 



SDFOO'^ 



Age [Myr] 
D.M. [mag] 
(Dist. [pc]) 



49.5±9.5 
7.72±0.10 
(350±16) 



43.0±16.0 
7.86±0.11 
(373J 



^-"-ID 



46.4 

7 7t;+oT04 
'■''^-0.05 

s) 



(355 



^ Intrinsic distance moduli were measured using an extinction of Ay = 0.54 
mag. 

^ Median of values from all membership catalogs for a given model and color. 
Uncertainties quoted are the range of possible values given the errors on the 
individual best-fit isochrone models. We did not include Ho isochrone mea- 
surement i n median and range calculation due to their large uncertainties. 
<= DAM97 -ID'Anto na fc M azzitellil l(T997i ) ; BCAH98 - IBaraffe et"an iflgM ) 
; SDFOO - ISiess et al. (2000J) 



models are i n agreement with the best-fit values for 
IBaraffe et alj model isochrones yielding systematically 
older ages (although they are within 1 ct of the other au- 
thors' models). For V — Ic, the derived age s for the three 
models are in close agreement, however, the lBaraffe et al.l 
models again yiel d slightly older ages. Furthermore, 
IBaraffe et alJ and iSiess et all models yield deriv ed DM 
that are larger than the ID'Antona fc Mazzitellil V — Ic 
models. 

Our results also show that all three models give sys- 
tematically older ages for the B — V CMDs compared to 
the V — Ic ones. The color-dependent age discrepancies 
seen are of order '--'10 Myr, and are on the same order of 
offsets seen in other young open clusters when comparing 
isochrone fits using two different colors (sec Naylor et ai] 
2002f). W e also note that for the IBaraffe et al.l and 



Siess et all models, there is a color-dependent offset of 
~0.2 magnitudes seen in the measured DM from the 
V -Ic CMDs. We find the ID'Antona fc MazzitelH mod- 
els provide the most consistent age and DM between the 
IC 4665 B-V andV-Ic CMDs. 

At this point, we caution the reader that the ob- 
served model- and color-dependent discrepancies in ages 
and DM are likely due to many different influences, 
and may not necessarily be representative of problems 
with low-mass stellar evolution models. Indeed, these 
three models include different physics that almost cer- 
tainly must lead to different predicted luminosities and 
effective temperatures at given ages (jSiess et al.l 120001 : 
iHillenbrand fc Whitell200^ for detailed discussion, see); 
however, the observed offsets may also be the result of 
systematics in the modeling process. For example, the 
ages derived from the B-V CMDs may be older due to 
the B-V color of late-type stars becoming nearly in- 
sensitive to changes in effective temperature, which con- 
sequentially makes it difficult to assign a PMS isochrone 
to such low-mass stars. In addition, the Pleiades tuning 
method may also introduce systematic errors into our 
analysis. For instance, Stauffer et al. (2007) find that 
low-mass stellar evolution models have trouble fitting the 



observed optical and infrared colors of low-mass members 
of the Pleiades. Because we calibrate our IC 4665 col- 
ors and magnitudes using an empirical fit to the Pleiades 
main-sequence, any such systematics in the Pleiades tun- 
ing method are included in our isochrone modeling proce- 
dure. More generally however, by employing the Pleiades 
tuning method we are assuming that the colors and mag- 
nitudes of the Pleiades main-sequence stars are represen- 
tative of the younger stars in IC 4665. Unfortunately, our 
IC4665 dataset is relatively small compared to other well 
studied open clusters and does not include a complete, 
contaminant- free PMS, which therefore may not allow 
us to disentangle all of the different color- and model- 
dependent systematics affecting PMS isochrones model- 
ing. 

Keeping this in mind, we would nevertheless still like to 
derive a single, best-fit PMS age an d distance for IC 4665 
for comparison purposes. Since the lJeffries et al.l catalog 
is most likely to contain the highest number of bona fide 
IC 4665 members, we use this catalog's resulting best-fit 
isochrone ages and DM as representative of the cluster. 
First, we calculate the inter-model error-weighted mean 
age and DM for the B-V and V — Ic colors, separately. 
The errors on these values are asymmetric, therefore the 
weights for our means are based on whether the age/DM 
lies above or below the straight average, e.g., if an age 
lies below the straight average, we then weight that data 
point using its upper error value. Using this approach, 
we find color-dependent ages of 30.6±3.5 and 49.3±9.0 
Myr, as well as intrinsic DM of 7.88±0.05 and 7.75±0.01 
mag (377±9 and 355±3 pc), for our B-V and V — Ic 
CMDs, respectively. Taking the weighted mean of these 
color-dependent values, we get a final PMS age and DM 
for IC 4665 of 35.8±9.3 Myr and 7.78±0.07 mag (360±12 
pc), for which the quoted errors are uncertainties in the 
mean. 

The entire photometric catalog of IC 4665 fields is plot- 
ted in Fig. \to\ in the form of V/V - Ic and V/B - V 
CMDs, where stars having photometric data consistent 
(within the errors) of being photometric members of 36 
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Fig. 10. — Intrinsic V — Ic and B — V color-magnitude diagrams for IC 4665 are plotted for our entir e photometric catalog. Stellar 
objects having colors and magnitudes consistent, within their photometric errors, of a 36 Myr and 360 pc'D'Antona & : Mazzitellil model 
isochrone (red solid line) are highlighted (blue diamonds), together with the locus of its associated equal- mass binary sequence (red dotted 
line). Photometric data have been dereddened using an Ay = 0.543, an E{B — V) = 0.174 and an E(V — Ic) = 0.2175. 



Myr and 360 pc iD'Antona fc Mazzitellil isochrones are 
highlighted. For each of these candidate photometric 
members of IC 4665, we list their astrometric and pho- 
tometric properties in Tabled This photometric mem- 
bership selection is based upon ID'Antona fc Mazzitellil 
because our model- and color-averaged age and DM are 
most closely matched to the values we get from the indi- 
vidual modeling of the V — Ic and B — V CMDs compared 
to the other theoretical PM S models gen erated by other 
authors (e.o.. lBaraffe et ahl fc lSiess et al.l ). Although this 
membership list most likely contains a high number of 
cluster members, its exclusive use in identifying mem- 
bers of IC 4665 is not recommended due to the non-zero 
probability of field star contamination in the catalog, es- 
pecially where the Galactic field star distribution inter- 
sects the cluster main-sequence close to {V — Ic)o ~ 1-00 
and {B - V)o ~ 1.25. 

Previously, we have shown that the CMDs of IC 4665 
members appear very similar to those of the ~35 Myr old 
NGC 2547 cluster (see Fig. [6]), which is now empirically 
supported by our PMS age determination for IC 4665. 
Moreover, the best-fit DM that we derive for the clus- 
ter is in excellent agreement with the trigonometric par- 
allax distance resulting from the recent r e-reduction of 
the H IP PARC OS astrometric dataset by Ivan LeeuwenI 
(|2009f ). which yields an intrinsic DM of 7.75±0.21 mag 
(355±35 parsecs). 

3.2. Upper-Main-Sequence Turn-off Age 

Comparing the positions of open cluster uppcr-main- 
sequences in color-magnitude space, .Mcrmilliod (,1981a,) 



determines that IC 4665 had a nuclear turn-off age of 
36 Myr. However, P93 suggest that the upper main- 
sequence for IC 4665 is best comparable to that of the 
~70 Myr open cluster a Persei, and if effects due to rapid 
rotation are taken into account, the upper main-sequence 
is most like the Pleiades CMD, suggestive of a cluster age 
near to -100 Myr. 

In light of this disagreement in nuclear turn-off age for 
IC 4665, we now use the r^ technique, outlined in Sec- 
tion 13.11 to measure an age and distance for only the 
high-mass cluster members. Because the main-sequence 
turn-off for IC 4665 is located at a V magnitude brighter 
than our CTIO photometric survey saturation limit, we 
utilize t he bright (F < 11 rnagnitu de) optical data from 
P93 and lMenzies fc Marangl (|199 6fl. Unfortunately, nei- 
ther catalog contained published individual uncertainties 
on their photometric data; therefore, we set the error in 
the B and V magnitudes to the level of observed scatter 
in P93 data in stars with V < 14 mag found in our cata- 
log comparison {AV = 0.07 mag and AB — V = 0.1 mag, 
see Section l2?3l) . Following [Navlor (2009) , we use upper - 
main-sequ ence mo dels from ILejeune fc Schaererl (|2001l ) 
and Girardi et al.l (|2002) , coupled with solar- metallicity 
bolometric corrections and effective temperature to color 
conv ersions, as wel l as co lor dependent reddening vectors 
from lBessell et all ()1998f ). We find upper-main-sequence 
ages for IC 4665 of 41±12 and 42±12 Myr and intrin- 
sic DM of 7.76±0.04 and 7.77±0.04 mag (356±10 and 
358±15 parsecs) for the aforementioned different high- 
mass stellar evolution models, respectively. We plot in 
Fig.HHthe photometric data used in our high-mass CMD 
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TABLE 6 
Catalog of IC 4665 Photometric Members 



Mosaic Object V 

Field Number 



B-V 



V-Ic 



Phot Mem^ 

[B-V] 



Phot Mem^ 
[V - Ic] 



Mem Key*" 



Fl 


427 


13.069±0.002 


0.816±0.005 


0.997±0.004 


Y 


F6 


347 


12.729±0.002 


0.835±0.004 


1.018±0.003 


N 


F6 


59 


14.475±0.006 


1.252±0.019 


1.274±0.009 


Y 


F6 


338 


12.095±0.001 


0.628±0.002 


0.754±0.002 


Y 



Y 
Y 
Y 
Y 



Li,RV,J, 

Li,XR,RP 

RV 



o 

> 



Note. — Only a few rows are displayed here to show the content of the table, the full table is available in 
electronic form. 

^ Point lies within its photometric uncertainty of the isochronc using the specified color. 

'' Ancillary membership criterion. See Section 13.11 for source references. Li = Lithium absorption, XR = 
X-ray emis sion, H = Ha emission, RV = Radial Velocity, RP = detected rapid rotation period, J = listed as 
member in ljeffries et al.l H2009al ) . and '...' = photometric member only. 
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Fig. 11. — High-mass (V<11 mag) color-magnitude diagr am us - 
ing photometric data from the optical catalogs of : Prossen 1)19931 ) 
and IMenzies &: Maraud H1996I V Also plotted are the best-fit 
isochrones (solid lines) from our upper-main-sequence T^-modeling 
using the two high-mass stellar models detailed in the legend, as 
well as their associated equal-mass binary sequences (dotted lines). 
The adopted uncertainty for the optical data is shown. 



modeling, as well as our two best-fit isochrones. We note 
that this age measurement is not thoroughly constrained 
due to the lack of high-mass cluster members tracing 
out IC 4665's upper-main-sequence; nevertheless, we find 
that the nuclear turn-off age for IC 4665 we have calcu- 
lated is similar to the previous estimate from IMermilliodI 
(|1981aD . and within the error budget, are in agreement 
with the a Persei-like age estimate from P93 (although 
quite discrepant from P93's Pleiades-like age estimate). 
iNavlor, (t2009,) recently investigated the differences ob- 



served in isochronal ages measured from upper- and 
lower-main-sequences of young (< 100 Myr) open clus- 
ters. They find a trend of PMS turn-on ages being ~1.5 
times smaller than nuclear turn-ojfa,ges for a selection of 
young open clusters. Interestingly, our current study on 
IC 4665 shows the upper-main-sequence isochrone mod- 
els have best-fit ages that are similar (difference of ~6 
Myr) to the PMS isochrone age while using the B — V 
low-mass CMD. However, if we use the PMS isochrone 
age from the V — Ic CMDs (30.6±3.5 Myr) or our inter- 
model F- 4 mean PMS age (35.8±9 Myr), we find -1.5 
times these PMS ages are formally within uncertainty 
of the nuclear turn-off age {i.e., 1.5 x 30.6 ± 3.5 Myr 
= 45.9 ± 3.5 Myr and 1.5 x 35.8 ± 9 Myr = 53.7 ± 9 Myr, 
compared to our turn-off age of 42 ± 12 Myr), corrobo- 
rating the lNavlod p009l findings. 

3.3. Lithium Depletion Boundary Age 

Using the lithium depletion boundary to date young 
open clusters provides a precise ('^ 10% uncertainty), 
nearly model-inse nsitive method to determine ages 
for o pen clusters () Burke et al.l 120041 : iJeffries &: Oliveiral 
|2005( ). It has been found that LDB ages for open 
clusters tend to be systematically older than upper- 
main-sequence turn-off age s (Stauffe r et al.|[l998l 119991 : 
iBarrado V Navascues et all 12004: .Burke et al.l 120041) : 
however, the discrepancy between LDB and turn-off a,ges 
has been found to be smaller for younger open clusters 
([Jeffries & Olivcira 2005; Navlor 2009) and may n ot exist 
for PMS isochrone ages ([Jeffries fc Oliveirall2005| ) . 

IC 4665 provides an interesting test for these age dis- 
crepancies in young open clusters because it is one of 
the few clusters with a LDB age in addition to well- 
defined PM S and upper-main-s equence populations. We 
recall that, iManzi et al.l (|2008f l measured a LDB age of 
28±5 Myr age for IC 4665, which is in agreement within 
the error budget, with our 35.8±9.3 Myr model-averaged 
PMS age determination. We also find the upper-main- 
sequence age for IC 4665 (42±12 Myr) is 1.5 times older 
than the LDB age, thereby contradicting the trend seen 
in other open clusters with LDB ages. We note, however, 
that the paucity of high-mass stars in PMS open clusters 
like IC 4665 reduces the fidelity of their turn-off ages, 
therefore this discrepancy in the observed turn-off-LDB 
age trend may turn out to be of little significance. 

3.4. X-ray Activity and Age 
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Fig. 12. — R atio of X-ray to bolometric luminosity plotted as a function of intrinsic V — Ic color for X-ray sources taken from 
IGiampapa et all {1998) with counterparts in our optical catalog for IC 4665 (blue asterisks). In each relevant panel, also plotted are 
X-ray data for NGC 2547 (red squares), IC 2391/IC 2602 (orange diamond s), a Pers ci (purple triangles), Blanco 1 (black crosscsV 
Pleiades (li ght blue triangl e s), and M 7 (green c ircles ). Data are taken fr om ljeffries et al. (2006), Stauffer et al. (1997), Randich ct ^ 
l|T996i). . Cargile et al.l 1)20091 ') . FStauffer et al.l H1994I '). and lProsser et al] 1)19951 ). respectively. Spectral-type ranges, as defined bv .BesscU et aP 
II1998I) . are plotted along the top axes. Intrinsic colors were calculated using E{B — V) values collated in the WEBDA database and an 
E{V - Ic) = 1.25 X E{B - V) relation. 



There is a very distinct, age- and mass-dependent mor- 
phology of X-ray activity observed in young open clusters 
that is fundamentally tied to the causal relationship be- 
tween age-dependent stellar rotation and the presence of 
a rotationally-induced, magneto-hydrodynamic dynamo 
in young stars (for review, see Giidel 2004). More specif- 
ically, the distribution of X-ray emission, typically ex- 
pressed as the distance-independent ratio of X-ray to 
bolometric luminosity (Lx/Ltoi), can be described by 
two regimes: the first, a regime where Lx/Lhoi more or 
less monotonically increases with decreasing stellar mass 
(and independently, with decreasing stellar period); and 
second, the so called saturated regime, where there is a 
plateau in X-ray emission at L^/Lhoi ~ 10"'^ regardless 
of ch anges in stellar mass or change in stellar ro tation 
rate (|Hempelmann et aLlfTOOl I James et"alll2000[ ). The 
physical cause of this saturated regime is stil l a matter of 
deba te (Charbonncau & MacGregor 19(3; 10 'dell et al.l 
Il995t [jardine fc Unruhl Il999( ): however, the transition 



point between these two regimes is empirically found 
to be a function of rotatio n period for a given stellar 
mass (|Pizzolato et al.l [20031 ) . and therefore, a function 
of age due to the time-dependence of stellar rot ation 
(|Barneall200.i [20071 : IMamaiek fc Hillenbrandll2008n . For 
example, in a young cluster like NGC 2547 (^^35 Myr), 
even the higher mass stars (^ G-type stars) are rotat- 
ing signific antly fast enough t o be in the saturated X- 
ray regime ([Jeffries et al.|[2000() , while in an older cluster 
like the Hyades (~700 Myr) the transition point between 
the two regimes i s found in the M dwarfs mass domain 
([Stern et al.lll995f ). Consequently, we can assign approx- 
imate ages of open clusters, in a relative manner, by com- 
paring their X-ray distributions to those of clusters with 
known ages, and comparing the observed color (i.e., stel- 
lar mass) at which there is a transition from the linearly 
increasing, color-dependent X-ray luminosity function to 
the saturated regime. 
In the ROSAT study of IC 4665 ([Giampapa et all 
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119980 , the authors suggest that IC 4665 has a X-ray dis- 
tribution indicative of a cluster which is older than the 
a Persei (~70 Myr), and very similar to the Pleiades 
(~120 Myr). They find that the color at which the 
IC 4665 X-ray distribution appears to transition to the 
saturated regime is redder (lower-mass) than a Persei, 
and is in agreement with the color of saturation for the 
Pleiades. Therefore, according to the age-activity rela- 
tionship, IC 4665 should have an age of ~120 Myr. This 
finding is formerly discordant with our current high- and 
low- mass isochronal analysis of IC 4665 (Sections [SH] and 
13. 2p . as well as its measured LDB age. 

In order to evaluate the validity of this X-ray age 
determination, we utilize our new optical photometry 
and sev eral X-ray survey s of other clusters published 
since the lGiampapa et al.l ([1998) study, including several 
with data from the new generation of X-ray telescopes 
(Chandra and XMM-Newton X-ray observatories). For 
ROSAT X-ray detected stars in IC 4665, we assign new 
colors (V — Ic) and calculate new bolometric luminosi- 
ties from our o ptical cat a log. We use the bolometric 
corrections from iJohnsonI ((1966!) for stars with intrinsic 

V — I n < 1.6 and the formalism given in iMonet et al.l 
(|1992l ) for stars with intrinsic V - Ic > 1.6. In Fig. [H 
we show our new Lx/Lboi values as a function of intrinsic 

V — Ic color for IC 4665, as well as the X-ray d istribution 
of the 35 Myr NGC 2 547 (IJeffries et al.l [2006l). 50 Myrs 
IC 2391/IC 2602 (.Stauffer et al.lll9970 . 70 Myr a Persei, 
100 Myr Blanc o llCargile et al.l l2009\ 120 Myr Pl eiades 
dStauffer erallll994[ ). and 250 Mvr M 7 (Prosser et al.l 

119951 ) open clusters . 

Looking at the overall morphology of X-ray emission 
in Fig. [121 we find that the X-ray distribution of IC 4665, 
and more specifically the color at which saturated stars 
are observed, is much more in agreement with the young 
NGC 2547 cluster than the older open clusters. Accord- 
ing to the relationship between age and activity, the sim- 
ilarity of the two cluster's X-ray activity distributions 
suggests that IC 4665 h as an age nearer to ~35 Myr 
([Jeffries fc Oliveirall2005[ ). thus is in agreement with our 
isochrone ages and the published LDB age. 

Although the IC 4665 data include our new SMARTS 
optical photometry, the onset of saturation in a Persei 
and the Pleiades is still observed to be very similar to 
IC 4665 (and furthermore NGC 2 547). In light of our 
larger cluster sample compared to iGiampapa et all we 
find that the X-ray emission distribution of the Pleiades 
and a Persei appears to be in disagreement with the 
overall observed act ivity-age morphology. We note that 
iCargile et al.l ([2009D similarly find the Pleiades X-ray dis- 
tribution is at odds with the approximately coeval open 
cluster Blanco 1, and speculate that this discrepancy 
is primarily due to scatter in the Pleiades photometric 
data, and is not attributable to the stellar properties of 
Blanco 1 members. 

3.5. Ha emission and Age 

Much like X-ray activity, studies of chromospheric ac- 
tivity have also been used as a good indicator of relative 
age for open cluste rs ([Prosser et al.lll99lLlStauffer et al.l 
[19911 : [Prossedll99"2l ). Through its causal relationship to 
the solar-type dynamo, the level of chromospheric emis- 
sion changes as a function of stellar rotation, and there- 
fore ultimately on age (,Mamaick fc Hillenbrand i2008i) . 



More specifically, as a star spins down with age, reduced 
emission is o bserved in ch romospheric lines, e.g., Ha, 
Ca II H & K ([Wilson|[1963|). 

Similar to the X-ray study of IGiampapa et al.l . P93 
claim that IC 4665 has an age of >100 Myr based on 
the upper envelope of its Ha equivalent width distribu- 
tion, as a function of pseudo-MK (pMK) spectral-type, 
being below that seen for the Pleiades, thus in direct 
conflict with our current isochrone dating analysis and 
its measured LDB age. P93 derived pMK spectral-types 
using a series of line ratios, mostly strong TiO features, 
measured from their low-dispersion spectra. Here, using 
our new photometric catalog and Ha data listed in Sec- 
tion 13. 1[ we reanalyze the Ha distribution in IC 4665 
to investigate the reported chromospheric activity-age 
estimate of '^100 Myr. In Fig. [131 left panel, we plot 
Ha equivalent width (EW Ha) taken from P93 as a func- 
tion of our new intrinsic V — Ic color. We also plot Ha 
data for open clusters NGC 2547 ( 35 Myr; I Jeffries et aLl 
l2000l) . IC 2391/IC 260 2 (50 Mvr: IStauffer et al.l ll997D 
and Blanco 1 (100 Mvr: lPanagi &: 0'delllll997l James et 
al. in prep). The level of Ha emission seen, as well as 
the color at which Ha is observed in emission, in IC 4665 
indicates that the cluster is younger than 50 Myr, and 
most likely has an age near (or younger) than that of 
NGC 2547 (35 Myr). 

When we plot t he Ha emission distribution of IC 4665 
and the Pleiades ([Stauffer fc Hartmann|[T987[ ) using V — 
Ic, instead of using the pMK spectral- type classification 
schema, as a proxy for effective temperature (Fig. [121 
right panel), we see that the two clusters have dis- 
similar morphologies. The onset of Ha in emission in 
IC 4665 is found at an earlier spectral-type compared 
to the Pleiades, thus suggesting that IC 4665 is younger 
than the ^120 Myr old Pleiades. This finding is con- 
trary to the chromospheric dating analysis in P93; how- 
ever, their Ha line measurements were performed over a 
much more limited range of spectral-types with the ear- 
liest being ~M0. With our more extensive sample of Ha 
measurements from three different sou rces (Prosser 1993; 
IMartin fc MontesI 119971 : IJeffries et al.ll2009aD . we clearly 
observe Ha emission in IC 4665 stars as early as spectral- 
type ~G0 ((V — Ic)o ^ 0.7), thereby suggesting a much 
earlier age estimation. 

4. SUMMARY & CONCLUDING REMARKS 

In this paper, we present the results of a new, homo- 
geneous, standardized BVIc photometric survey of the 
central region of the young open cluster IC 4665. This 
dataset improves upon the previous work of P93, by in- 
corporating a high-fidelity standardization procedure in- 
volving observations of multiple standards taken contem- 
poraneously with our IC 4665 data. Our full photometric 
catalog contains 7359 detected point sources, of which 
all but 452 are matched within 3.0 arcseconds of 2MASS 
sources. The internal stability of our photometric sys- 
tem, as determined by a nightly control field, is 3 — 4% 
for objects with B, V, and Ic < 16 magnitudes. Using 
our best-fit age and distance from our isochrone model- 
ing (see l3.1.2[) we have identified 382 IC 4665 photometric 
cluster member from our photometry; 324 and 176 from 
the V—Ic and B—V CMDs, respectively. We provide this 
extensive membership catalog due to it likely containing 
a large number of cluster members, we advise that it not 
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Fig. 13. — Equivalent width (EW) of the Balmer series Ha spectral line is plotted versus intrinsic V — Ic color for matched optical 
counterparts in our IC 4665 catalog (blue asterisks), as well as data for (left panel) NGC 2547 (red squares), IC 2391/IC 2602 (orange 
diamonds), B lanco 1 (black crosses), and (right panel) the Ple i ades. Hg data are tak en fr om Prosscr (1993); Ma rtin fc Montes (1997); 
IJeftries et all (|2009al ). I Jeffries ct al. (2000), Stauffcr et al. ( 19971. Panagi fc 0'delll ((T997l ). and [S'tauffer fc Hartmarml ||T987D for each cluster, 
respectively. Spectral-type ranges, as defined by BesscU ct al. (1998), arc plotted along the top axes. We deredden the V — Ic colors using 
the cluster reddening coefficients listed in WEBDA, and defined positive EWua to be in emission. 



be used exclusive to identify clusters members due to a 
high level of field star contamination in IC 4665. 

In order to refine our photometric membership catalog 
into a high-fidelity membership list for IC 4665, we cor- 
relate the source positions in our new BVIc photometric 
survey with those in various extant membership cata- 
logs. This includes membership based on kinematic prop- 
erties, lithium absorption, rapid rotation, and chromo- 
spheric/coronal emission. We employ these membership 
catalogs with our updated photometry to reevaluate the 
age of IC 4665 using various dating techniques. New pre- 
main-sequence turn-on and upper-main-sequence turn- 
off ages were calculated using r^ modeling, a least- 
squares fitting procedure well suited for modeling color- 
magnitude diagrams due to it taking into account corre- 
lated errors in age and distance. We find a best-fit PMS 
isochrone age and distance of 35.8±9 Myr and 360±12 
pc, and a turn-off age and distance of 42±12 Myr and 
357±13 pc, where these values are weighted averages and 
the errors are uncertainties on the means. 

Definite model- and color-dependent offsets are seen 
over a range of ^^10-20 Myr in age and ~20 parsecs 
in distance, particularly in the results from the PMS 
isochrones. This large spread in ages and distances, as 
well as the systematic offsets seen between colors and 
models, are likely due to limitations in the individual 
membership catalogs {e.g., a narrow range of colors for 
a given membership catalog produces large uncertain- 
ties in its derived age and distance) and/or systematic 
differences in the predicted stellar parameters from the 
different sets of models. At this point, how these sys- 
tematic model- and color-dependent offsets infiuence the 
measured ages and distances cannot be completely de- 
ciphered primarily due to a limited number of known 



clusters members and a relatively high non-member con- 
tamination level of IC 4665. 

In addition, we also compared the distribution of mag- 
netic activity in IC 4665 stars to other well studied open 
clusters. The coronal X-ray and chromospheric Ha emis- 
sion observed in IC 4665 is most similar to 30-40 Myr old 
open clusters {e.g., NGC 2547), thereby providing fur- 
ther confirmation on our measured turn-on and turn-off 
ages. This result also provides further evidence against 
previous claims that IC 4665 has an observed activity 
distribution similar to the Pleiades. Previous IC 4665 
activity-age studies have been reliant on relative distri- 
bution comparisons with clusters that we find do not fit 
into the observed global trends seen in larger samples of 
open clusters, particularly the X-ray emission distribu- 
tion of a Persei and the Pleiades. 

While it is clear that we do see slight offsets and dif- 
ferences in the various dating techniques we have em- 
ployed in this paper, the majority of our measured ages 
for IC 4665 are consistently between ~30-45 Myr. Inter- 
estingly, the lack of large systematic differences (factors 
of >1.5 in age) between dating techniques is not typi- 
cal for open clusters; for exa mple, published ages for the 
Pleia des range from 70 Myr ()Mermilliodlll981b[ ) to 150 
Myr (jBurke et al.ll2004D . Understanding why IC 4665 
does not suffer from these systematic offsets should allow 
for a critical comparison between the different methods 
available to determine stellar ages. 
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